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We study the effect of tensor correlations on single-particle and collective 
states within Skyrme Hartree-Fock and RPA model. Firstly, We study the 
role of tensor interactions in Skyrme effective interaction on the spin-orbit 
splittings of N=82 isotones and Z=50 isotope. The isospin dependence of 
the shell structure is well described as the results of the tensor interactions 
without destroying good properties of the binding energy and the rms charge 
radii of the heavy nuclei. Secondly, We performed self-consistent HF+RPA 
calculations for charge exchange 1 + states in 90 Zr and 208 Pb to elucidate the 
role of tensor interactions on spin dependent excitations. It is pointed out 
(N| ■ that Gamow-Teller(GT) states can couple strongly with the spin-quadrupole 

(SQ) 1 + states in the high energy region above E x =30 MeV due to the tensor 
interactions. As the result of this coupling, more than 10% of the GT strength 
is shifted to the energy region above 30 MeV, and the main GT peak is moved 
2 MeV downward. 
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1. Introduction 

One of the current topics in nuclear physics is the role of the tensor interactions 
on the shell evolution of single-particle states and on spin-isospin excitations. The 
importance of tensor interactions on nuclear many-body problems was recognized 
more that 50 years ago. Especially, it plays the essential role to make the binding 
f-^ ; systems such as the deuteron. As is shown Fig. [U the tensor interaction acts on the 

spin triplet state (S=l) of two nucleon system. When a proton and a neutron are 
. aligned in the direction of spins, the deuteron gets an extra binding energy by the 

tensor interaction, 

V T = f(r)S 12 (1) 
On . 

since f(r) is negative and Su = 3(<?i • f)(a 2 • r) — ff\ • a 2 =2. On the other hand, if 
a proton and a neutron are perpendicular to the spin direction, the deuteron will 
loose the binding energy since £12 = — 1 and cannot make a bound state. Thus, the 
deuteron becomes a strongly deformed prolate shape and makes a binding system 
only when the tensor correlations are taken into account. 

After these findings, the importance of tensor interactions has been widely recog- 
nized in nuclear many-body systems, especially in light nuclei. The role of the tensor 
interactions in the Hartree -Fock calculations was firstly discussed by Stancu et al., 
thirty years ago [1]. However serious attempts have never been performed until 
very recently [2-6]. The importance of the tensor correlations on the mean field 
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Figure 1: Tensor interaction. See the text for details. 



revived by the study of the shell evolution of heavy exotic nuclei [7]. In this paper, 
I will summarize formulas of Skyrme tensor interactions in Section 2. The isotope 
dependence of the shell structure of Z=50 isotopes and N=82 isotones are studied 
by using the Hartree-Fock (HF)+BCS calculations in Section 3. The role of tensor 
interactions on spin-isospin excitations will be examined in Section 4. Summary is 
given in Section 5. 



2. Skyrme tensor interaction 

The tensor force was considered in the Skyrme-Landau parameterization and the 
sum rules of electromagnetic transitions in Ref. [8]. However, the tensor force was 
essentially dropped in most Skyrme parameter sets which have been used widely in 
nuclear structure calculations. Recently, in Ref. [3], a Skyrme interaction was fitted 
including the tensor contribution. Then, tensor terms were added perturbatively in 
Refs. [2] and [4] to the existing standard parameterizations SIII [9] and SLy5 [10], 
respectively. Eventually, several new parameter sets have been fitted in Ref. [5] 
and used for systematic investigations within the Hartree-Fock-Bogoliubov (HFB) 
framework. The inclusion of tensor terms in the Skyrme HF calculations achieved 
considerable success in explaining some features of the evolution of single-particle 
states [2,6]. The Skyrme tensor interaction is given by the triplet-even and triplet- 



odd tensor zero-range tensor parts, 

v T = \{{{<J l -k'){a 2 -k')-\{a 1 -a 2 )k' 2 ]8{r l -r 2 ) 
+ 5(n - r 2 )[(a 1 • k)(a 2 ■ k) - ^ • a 2 )k 2 ]} 

+ U{{a x ■ k')8(n - r 2 )( ( r 1 • k) - ~ (<7i • <7 2 )[fc' • <5(r 2 - r 2 )A:]} (2) 

where the operator k = (Vi — V2)/2z acts on the right and k' = — (Vi — V2)/2i 
on the left. The coupling constants T and [/ denote the strength of the triplet-even 
and triplet-odd tensor interactions, respectively. We treat these coupling constants 
as free parameters in the following study. The tensor interactions ([2D give the 
contributions to the binding energy and the spin-orbit splitting proportional to the 
spin density 



jiUi + 1) - k{k + 1) - \ 



R-{r) (3) 



8E = -a(Jl + J 2 )+pJ n J p . (4) 



where i = n, l,j runs over all states and q = 0(1) is the isospin quantum number for 
neutrons (protons). The vf is the occupation probability of each orbit determined 
by the BCS approximation and Ri(r) is the HF single-particle wave function. It 
should be noticed that the exchange part of the central Skyrme interaction gives the 
same kind of contributions to the total energy density. The central exchange and 
tensor contributions give the extra terms to the energy density as 

1 

2 ( 

The spin-orbit potential is then expressed to be 

s -°- 2r \ dr dr J \ r r J 

where the first term on the r.h.s comes from the Skyrme spin-orbit interaction and 
the second term include both the central exchange the tensor contributions a = 
a c + (%t and /3 = /3 C + Pt- I n Eq. g = 0(l)is assigned for neutrons (protons). 
The central exchange contributions are given by 

!/ x 1/ 

OLC = o (*1 ~ *2) - o \t\Xl + t 2 X 2 ) 
o o 

/?c = (6) 

o 

where the parameters are defined in ref. [11]. The tensor contribution are expressed 
as 

5 

a T = -U 

Pt = Y^T + U). (7) 



Effect of tensor interaction on spin-orbit splitting 
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Figure 2: Tensor and Spin-orbit interaction. See the text for details. 

Before going to detailed study, let us mention two important features of the tensor 
and the central exchange contributions in Eq. ([5]) to the spin-orbit splitting. First 
point is that the mass number dependence of the the first and second terms in Eq. (ED . 
Since the Skyrme spin-orbit force Wo gives the spin-orbit splitting proportional to 
the derivatives of the densities, the mass number dependence is very modulate in 
heavy nuclei. On the other hand, the second term in Eq. (|5j) depends on the spin 
density J q which has essentially no contribution for the Is closed shell. The spin 
density will increase exactly proportional to the number of particles in the open 
shells if one of the spin-orbit partner is only active. Moreover, the sign of the J q 
will change depending upon which orbits are involved in the active shell orbits, i.e., 
the orbit j> — I + 1/2 gives a positive J q value while the orbit j< = / — 1/2 gives 
a negative J q . This means that the spin-orbit energy will change in the opposite 
direction according to which orbit is occupied in the open shell nuclei. 
We fit the two parameters T and U (equivalent ly ctr an d fir) using the recent 
experimental data of TV =82 isotones and Z = 50 isotopes. We keep the central 
part of Skyrme interaction as that of SLy5. The central exchange interactions are 
a c =80.2MeV-fm 5 and fi c = -48.9 MeV-fm 5 for SLy5. The optimal parameters a T 
and fir are determined to be (a^, fir) = (— 170,100)MeV-fm 5 . We examine detailed 
properties of the tensor interactions by using our parameter sets. In Fig. El we 
consider a nucleus where the last occupied orbit is the neutron j> = 1 + 1/2 (for 
example, imagine 90 Zr or 48 Ca). In. Eq. ([5]), the ar increase the spin-orbit splitting 
because of positive J q=0 of j> = I + 1/2 orbit, while fir has the opposite sigh and 
decrease the splitting for positive J q =o. These effects are illustrated in Fig. [2] where 
the neutron spin-orbit splitting is increased, while the proton one is decreased. 
In Fig. [3l the energy differences for proton single-particle states Ae(hu/2 — c/7/2) 
of Z=50 isotones are shown as a function of the neutron excess (TV — Z). The 



original SLy5 interaction fails to reproduce the experimental trend qualitatively and 
quantitatively. Firstly, the energy differences of the HF results are much larger than 
the empirical data. Secondly, the experimental data decrease, the neutron excess 
decreases and reach about 0.5MeV at the minimum value. On the other hand, the 
energy differences of the original SLy5 increase as the neutron excess decrease and 
has the maximum at around (N — Z)—20. We studied also several other Skyrme 
parameter sets and found almost the same trends as those of SLy5. 
The tensor central exchange interactions are included in the results marked by open 
circles in Fig. [3l We can see a substantial improvement by introducing the tensor 
interactions. The set (a^, Pt) = (— 170,100)MeV-fm 5 gives a fine agreement with 
the empirical data from (TV — Z) = 20 ~ 32 quantitatively and qualitatively. 




Figure 3: Comparison of energy difference between pairs of single-particle lg7/2 and 
lhn/2 proton states of Z=50 isotones. The calculations are performed with and 
without tensor and central exchange terms in the spin-orbit potential (|5j) top of 
SLy5 parameter set. The experimental data are taken from ref. [12]. See the text 
for details. 

The HF+tensor results can be qualitatively understood by general argument as 
follows. Firstly, the strength a change the depth of the proton spin orbit potential. 
For Z=50 core, only the proton g 9 / 2 proton orbit dominates the spin density J p in 
Eq. ((3j) so that with the negative c^r= value the spin-orbit splittings are increased. 
Thus, the g 9 / 2 protons increase the spin-orbit splitting between proton (g 9 / 2 — #7/2) 
orbits and that of proton /113/2 — ^11/2 orbits by qit effect. As a net effect, the 



Ae(/&n/2 _ 97/2) protons decreases substantially. 

Next let us study the (N-Z) dependence where (3t = plays the essential role. In the 
case of Ae(/in/ 2 — 97/2) protons on Z=50 core from N-Z=(6-14), the g 7 / 2 neutron 
orbit is gradually filled. Then the /3r = 100MeV-fm 5 gives a negative contribution 
to the spin-orbit potential ([5]) and increases the spin-orbit splitting. Therefore the 
energy difference Ae(/in/ 2 — 97/2) is further decreasing. From N— Z=(14 to 20), the 
S1/2 and ds/2 neutron orbits are occupied. In this region the spin density is not so 
much changed since Si/ 2 has zero contribution. For N-Z= (20-32), the hi 1/2 orbit 
is gradually filled. Then this orbit gives a positive contribution to the spin-orbit 
potential ([5]) ,i.e., the the spin-orbit splitting is decreasing. Then the Ae turns 
out to be increasing. The magnitude of f3 term determines the slope of the (N-Z) 
dependence so that a larger (3 gives a steep slope. 
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Figure 4: Comparison of energy difference between pairs of single-particle H13/2 and 
lhg/2 neutron states of N=82 isotopes. The calculations are performed with and 
without tensor and central exchange terms in the spin-orbit potential ([5]) top of 
SLy5 parameter set. The experimental data are taken from ref. [12]. See the text 
for details. 

For the results of Sb-isotones in Fig. [H the A e(/in/ 2 — #7/2) f° r N=82 core is 
plotted as a function of neutron excess. Essentially, the same argument can be 
applied for the neutron energy difference A e(ii 3 / 2 — ^9/2) for N=82 core as that for 
the A e(/iii/ 2 — #7/2) for Z=50 core. The last occupied neutrons in /in/2 increase 
the neutron spin-orbit splitting so that the A e(/in/ 2 — #7/2) become substantially 



smaller by the qit effect in Eq. (ED- The isotope dependence is again explained by 
the ft? effect. The lg7/2 and 2d 5 / 2 are almost degenerate above the last occupied 
proton orbit lgg/2 of Z=50 core. The two proton orbits, d 5 / 2 and g 7 / 2 , have opposite 
effects on the spin orbit potential (|5j). The occupation probability is larger for larger 
j orbit so that lg7/2 plays more important role on the spin orbit potential due to 
the tensor interaction in nuclei with (N-Z)=(32-18) for N=82 isotones. Namely the 
neutron spin orbit splitting is larger for these isotones so that the z 13 / 2 orbit is down 
and the h 9 / 2 is up. These changes make the energy gap A e(ii 3 / 2 — h 9 / 2 ) smaller for 
the nuclei from (N-Z)= 32( 132 Sn) to (N-Z)=18( 146 Gd). Thus the role of the triple- 
even and triplet-odd tensor interactions are clearly shown in Figs. [3] and H] for both 
proton and neutron spin-orbit splittings. 

The role of the tensor interaction due to the (3 term is essentially expected from the 
discussion by Blatt-Weisskopf [13] for the deuteron. The role of is new and has 
not been examined in a quantitative way in the mean field calculations since this term 
comes from the triplet-odd tensor interaction. The triplet-odd tensor interaction was 
not included in the studies of refs. [7,13]. Recently, Brown et al. studied the Skyrme- 
type tensor interactions in 132 Sn and 114 Sn based on the parameter set SkX. They 
took both the positive and negative values in the HF calculations and concluded 
that the negative ax value gives a better agreement with the experimental data. 
This is consistent with the present systematic study of Z=50 isotopes and N=82 
isotones with HF+BCS model. 

3. Tensor effect on Gamow- Teller states 

There has been no RPA or QRPA (Quasiparticle Random Phase Approximation) 
program available to study the effect of the tensor terms on the excited states of 
nuclei until very recently. The tensor terms of the Skyrme effective interaction was 
firstly introduced in the self-consistent HF plus RPA calculations [14], in particular, 
in the GT transitions, which should be affected because of the fact that the corre- 
sponding operator is spin-dependent. In the study of GT transitions, the quenching 
problem is of some relevance. The experimentally observed strength from 10 to 20 
MeV excitation energy (with respect to the ground state of the target nuclei) is 
about 50% of the model-independent non-energy weighted sum rule (NEWSR) [15]. 
It would be very interesting to study whether the tensor force has an effect in shift- 
ing the strength already at one particle-one hole (lp-lh) level. Coupling the GT 
with two particle-two hole states is essential to describe the resonance width but it 
is not expected to affect strongly the position of the main GT peak; the effect of the 
tensor force in connection with the 2p-2h coupling was studied in Ref [16]. 
It should be noted that J q gives essentially no contribution in the spin-saturated 
cases. Therefore, we choose 90 Zr and 208 Pb as examples to be calculated. 90 Zr is a 
proton spin-saturated nucleus, with a spin-unsaturated neutron orbit lgg/2- 208 Pb 
is chosen as it is not saturated either in protons or neutrons. The two examples 



should allow elucidating separately the role of triplet-even and triplet-odd terms. 
Since the tensor force is spin-dependent and affects the spin-orbit splitting, the spin 
mode is very likely to receive strong influence. The operator for GT transitions is 
defined as 

o GT± = J2 f ±< ( 8 ) 

im 

in terms of the standard isospin operators, t± = ^(t x ±it y ). In the charge-exchange 
RPA, the t- and t+ channels are coupled and the corresponding eigenstates emerge 
from a single diagonalization of the RPA matrix. 
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Figure 5: The GT_ strength in 90 Zr and 208 Pb. The RPA results are displayed, by 
smoothing them with Lorentzian function having 1 MeV width. As explained in the 
text, result labeled by 00 corresponds to neglecting the tensor terms in both HF 
and RPA; 10 corresponds to including the tensor terms in HF but neglecting them 
in RPA; finally, 11 corresponds to including the tensor terms in both HF and RPA. 
The arrow denotes the experimental energy. See the text for details. 

The GT_ strength distributions in 90 Zr and 208 Pb are shown in Fig. El The cal- 
culated results are smoothed by averaging the sharp RPA peaks with Lorentzian 
function weighting function having 1 MeV width. The tensor force affects these 
results in two ways. Firstly, it changes the single-particle energies (s.p.e.) in the 
HF calculation; secondly, it contributes to the RPA residual force. We do three 
different kind of calculations to analyze separately these effects. In the first one, the 
tensor terms are not included at all. In the second one, we include tensor terms in 
HF but drop them in RPA. This calculation is not self-consistent, but it displays 
the effects of changes in single-particle energies on the strength distribution. In 
the last one, the tensor terms are included both in HF and RPA calculations. For 
simplicity, results of the three categories of calculations are labeled by 00, 10 and 
11, respectively. 



Table 1: Values of the NEWSR m_(0) and EWSRs m_(l) for 90 Zr and 208 Pb in 
different excitation energy regions. The two-body spin-orbit interaction is included 
in HF but neglected in RPA calculation. The results labeled by 00 correspond to 
neglecting the tensor terms both in HF and RPA; 10 corresponds to including the 
tensor terms in HF but neglecting them in RPA; 11 corresponds to including the 
tensor terms both in HF and RPA. See the text for a discussion of the effects of the 
tensor terms. 
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2212.8 


18.8 
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10 


127.38 


3.68 


2176 


93 


2269 


21 




11 


114.10 


16.58 


1658 


694 


2370 


19.3 



We have evaluated the amounts of NEWSR m_(0) and EWSR m_(l) in different 
excitation energy regions, and listed them in Table [TJ The EWSR in the energy 
region below 30 MeV (where the one particle-one hole transitions are located) is 
decreased, after the inclusion of the tensor term. From Table [U we also see that an 
appreciable amount of EWSR is shifted from the lower energy region (0-30 MeV) 
to the higher energy region (30-60 MeV) by including tensor terms in HF plus RPA 
calculations. 

We also calculated the values of NEWSR in the 0-30 MeV and 30-60 MeV energy 
regions for 90 Zr and 208 Pb. When the tensor is not included in the residual interaction 
(i.e., the calculations labeled by 00 and 10), the values of NEWSR in the energy 
region between 30-60 MeV for both 90 Zr and 208 Pb are small only few percent of the 
NEWSR(Fig, El But in the case 11, about 10% of NEWSR is shifted from the lower 
energy region to the higher energy region (Corresponding 25% and 29% of EWSR in 
90 Zr and 208 Pb, respectively). Moreover, we can see that essentially no unperturbed 
strength appears in this region (see the Fig. [5|). This means that including tensor 
terms in simple RPA calculation shifts about 10% of the GT strength to the energy 
region 30-60 MeV. While 2p-2h couplings will increase further these high energy 
strength, we would like to stress that the tensor correlations move substantial GT 
strength from the low energy region 0-30 MeV to the high energy region 30-60 MeV 
even within the lp-lh model space. 

In 90 Zr, one can notice that the GT strength is concentrated in two peaks in the 
region below 30 MeV. There are only two important configuration involved which 
are (nlgg/2 — vlg^ 2 ) and (7iTg 7 /2 — vlg^J^ (see the left panel of Fig. [5]). When the 
tensor term is included only in HF and neglected in RPA, the centroid in the energy 
region of 0-30 MeV are moved upwards by about 1.5 MeV, and the high energy peak 



at Ex ~ lQMeV is moved upwards by only 0.5 MeV, as compared with the results 
without tensor term. When the tensor term is included both in HF and RPA, the 
centroid of the GT strength in the energy region 0-30 MeV is moved downwards 
by about 1 MeV, and the high energy peak is moved downwards about 2 MeV, as 
compared with the results obtained without tensor term. Including tensor terms 
in RPA makes the two main separated peaks closer (this situation also happens 
for 48 Ca). This result can be attributed from the HF and RPA correlations of the 
tensor term. When the vlg§/2 orbit is filled by neutrons, the tensor correlations give 
a quenching on the spin— orbit splitting between nlg$/2 and 7rlg7/ 2 orbits so that 
the unperturbed energies of the two main p — h configurations (7iTg7/ 2 — vlg^J^) 
and (7iTg 9 /2 — vlg^j^) are closer in energy as is shown in Fig. El The RPA results 
in Fig. [5] with labeled by (00) and (10) reflect these changes of HF single particle 
energies due the tensor correlations and the energy difference between two peaks is 
narrower. 

In 208 Pb, from the right panel of Fig [5] we see that the GT strength is concentrated 
in two peaks in the low energy region of 0-30MeV for all 00, 10 and 11. There 
are eleven important configurations which do contribute to these peaks. When the 
tensor terms are only included in HF and neglected in RPA, the centroid of these 
peaks is moved upwards about 0.5 MeV, and the higher energy peak at E x ^18MeV 
is also raised by about 0.8 MeV. When the tensor terms are included in both HF and 
RPA calculation, the centroid of these peak moves downward by about 1.5 MeV, and 
the higher energy peak moves also downwards by about 3.3 MeV, compared with the 
result obtained without tensor terms. By including tensor terms in RPA calculation, 
the GT strengths in the energy region of 30-60 MeV are increased substantially by 
the shift of the strength in the energy region of 0-30 MeV through the tensor force. 

4. Summary 

We study the effect of tensor correlations on single-particle and collective states 
within Skyrme Hartree-Fock and RPA model. Firstly, We study the role of tensor 
interactions in Skyrme effective interaction on the isospin dependence of spin-orbit 
splittings in N=82 isotones and Z=50 isotope. The different role of the triplet-even 
and triplet-odd tensor forces is elucidated by analyzing the spin-orbit splittings 
in these nuclei. The experimental isospin dependence of these splittings cannot be 
described by HF calculations with standard Skyrme forces, but is very well accounted 
for when the tensor forces are introduced, the GT excitations in 90 Zr and 208 Pb in 
the HF plus RPA framework with a Skyrme interaction SIII. If the tensor term is 
included in both HF and RPA, the centroid of G-T strength in the energy region 
below 30 MeV is moved downwards by about 1 MeV for 90 Zr and 3.3 MeV for 208 Pb. 
At the same time, the dominant peak at E x ^16MeV(18MeV) in 90 Zr( 208 Pb) is also 
moved downwards by about 2 MeV(3MeV). It is pointed out for the first time that 



about 10% of NEWSR is moved in the high energy region of 30-60 MeV by the 
tensor correlations in RPA even within lp — lh model space. It was pointed out 
recently that the high energy GT strength is shifted by the coupling between GT ana 
spin-quadrupole states due to the tensor correlations which has the intrinsic strong 
coupling [14] It is interesting to point out that the main GT peak, contrarily, gets 
the energy shift downward because of the peculiar feature of the tensor correlations. 
The tensor interaction is spin-dependent, so we expect that it can have important 
effects not only on the GT transitions, but on spin-dipole and other spin dependent 
excitation modes as well. 
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